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Abstract 

It is important to have effective methods for simplifying 3-manifold triangulations with- 
out losing any topological information. In theory this is difficult: we might need to make 
a triangulation super-exponentially more complex before we can make it smaller than its 
original size. Here we present experimental work that suggests the reality is far different: 
for an exhaustive census of 81 800 394 one- vertex triangulations that span 1 901 distinct 
closed orientable 3-manifolds, we never need to add more than two extra tetrahedra, we 
never need more than a handful of Pachner moves (or bistellar flips), and the average num- 
ber of Pachner moves decreases as the number of tetrahedra grows. If they generalise, these 
extremely surprising results would have significant implications for decision algorithms and 
the study of triangulations in 3-manifold topology. 

Key techniques include polynomial-time computable signatures that identify triangu- 
lations up to isomorphism, the isomorph-free generation of non-minimal triangulations, 
theoretical operations to reduce sequences of Pachner moves, and parallel algorithms for 
studying finite level sets in the infinite Pachner graph. 

ACM classification F.2.2; G.2.1; G.2.2; D.1.3. 

Keywords Triangulations, 3-manifolds, Pachner moves, isomorphism signatures, iso- 
morph-free enumeration, 3-sphere recognition 



Journal version 

This is the journal version of 9 , which was presented at the 27th Annual Symposium 
on Computational Geometry. This journal version contains significant new material. 

The study has been expanded from just 3-sphere triangulations to include all closed 
prime orientable 3-manifolds, it examines average-case as well as worst-case behaviour, and 
it also analyses moves that connect distinct minimal triangulations of the same 3-manifold. 
Section |4] contains several new algorithms, through which the loose upper bounds of [9] 
have been replaced with smaller bounds, many of which are now tight. The analysis of 
pathological cases in Section [5] and the detailed specification of isomorphism signatures in 
the appendix are both new. The final discussion in Section |S] is significantly richer, and 
raises new issues involving generic complexity. 
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1 Introduction 



Triangulations of 3-manifolds are ubiquitous in computational knot theory and low-dimensional 
topology. They are easily obtained and offer a natural setting for many important algorithms. 

Computational topologists typically allow triangulations in which the constituent tetrahedra 
may be "bent" or "twisted" , and where distinct edges or vertices of the same tetrahedron 
may even be joined together. Such triangulations (sometimes called generalised triangulations) 
can describe rich topological structures using remarkably few tetrahedra. For example, the 3- 
dimensional sphere can be built from just one tetrahedron, and more complex spaces such as 
non-trivial Seifert fibred spaces can be built from as few as two |25| . 

An important class of triangulations is the one-vertex triangulations, in which all vertices of 
all tetrahedra are identified together as a single point. These are simple to obtain [HI [24], and 
they are often easier to deal with both theoretically and computationally [6l \17\ [24] . 

Keeping the number of tetrahedra small is crucial in computational topology, since many 
important algorithms are exponential (or even super-exponential) in the number of tetrahedra 
[5l [6]. To this end, topologists have developed a rich suite of local moves that allow us to 
change a triangulation without losing any topological information [2l I26j . The ultimate aim is 
to simplify the triangulation, i.e., reduce the number of tetrahedra, although the triangulation 
might (temporarily) need to become more complex along the way. 

The most basic of these moves are the four Pachner moves (also known as bistellar moves). 
These include the 3-2 move (which reduces the number of tetrahedra but preserves the number 
of vertices), the 4-1 move (which reduces both numbers), and also their inverses, the 2-3 and 
1-4 moves. It is known that any two triangulations of the same closed 3-manifold are related by 
a sequence of Pachner moves ^34j- Moreover, if both are one- vertex triangulations then in most 
cases we can relate them using 2-3 and 3-2 moves alone [24] . 

However, little is known about how difficult it is to simplify a triangulation, or to convert 
one triangulation into another, using Pachner moves. In a series of papers, Mijatovic develops 
upper bounds on the number of moves required for various classes of 3-manifolds [29[ I30[ I3H 132] . 
All of these bounds are super-exponential in the number of tetrahedra, and some even involve 
exponential towers of exponential functions. For simplifying one-vertex triangulations using 
only 2-3 and 3-2 moves, no explicit bounds are known at all. 

Simplification is tightly linked to the important recognition problem, where we are given 
an input triangulation T and a target 3-manifold M., and asked whether T triangulates M.. 
The recognition problem is decidable but extremely difficult. A general algorithm comes as a 
consequence of Perelman's celebrated proof of the geometrisation conjecture [TT, but due to its 
intricate and multi-faceted nature, the algorithm remains computationally intractable and no 
explicit bound on its running time is known. 

Some special cases of the recognition problem are more approachable. A notable case is 
3-sphere recognition (where M = S^): this plays an key role in other important topological 
algorithms such as connected sum decomposition [18l[20] and unknot recognition [15], and is also 
important for computational J^-manifold topology. The original 3-sphere recognition algorithm 
of Rubinstein ^35, has been improved significantly over time [6. 7, ,18, , 36, , and although it remains 
worst-case exponential, it is now highly effective in practice for moderate-sized problems [7]. 

We can use Pachner moves to solve recognition problems in two ways: 

• For the classes of manifolds M studied by Mijatovic (29l [30] [31] [32], and in particular 
the 3-sphere |29j . Pachner moves give a direct recognition algorithm: select a well-known 
"canonical" triangulation 7^ of , and try to convert the input triangulation T into Tc by 
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testing every possible sequence of Pachner moves up to Mijatovic's upper bound. Return 
"true" if and only if a successful conversion was found. 



• For all manifolds M, Pachner moves also give a hybrid recognition algorithm: begin with 
a fast and/or greedy procedure to simplify T as far as possible within a limited number 
of moves. If we reach a well-known canonical triangulation of then return "true"; 
otherwise run the full recognition algorithm (such as Rubinstein's algorithm for the case 
A4 = S^) on our new, and hopefully simpler, triangulation. 

Direct algorithms are, at present, completely infeasible: Mijatovic's bounds are super- 
exponential in the number of tetrahedra, and the running times are super-exponential in Mija- 
tovic's bounds. Even for the trivial case of 3-sphere recognition with one tetrahedron, a direct 
algorithm must test all possible sequences of ^ 2.4 x 10^°^'' Pachner moves. 

The hybrid method, on the other hand, is found to be extremely effective in practice. Expe- 
rience with 3-sphere recognition software [3] suggests that when T is indeed the 3-sphere, the 
greedy simplification almost always gives a canonical triangulation, which means that the slower 
Rubinstein method is almost never required. 

In the context of Mijatovic's results, this effectiveness of the hybrid method is unexpected, 
and forms a key motivation for this paper. More broadly, the aims of this paper are: 

(i) to measure how difficult it is in practice to relate two triangulations of a 3-manifold using 
Pachner moves, or to simplify a 3-manifold triangulation to use fewer tetrahedra; 

(ii) to understand why greedy simplification techniques work so well in practice, despite the 
prohibitive theoretical bounds of Mijatovic; 

(iii) to investigate the possibility that Pachner moves could be used as the basis for a direct 
3-sphere recognition algorithm that runs in sub- exponential time. 

We restrict our attention to closed prime orientable 3-manifolds, as well as the important 
case A4 = (which is not prime). We also restrict our attention to one- vertex triangulations 
with 2-3 and 3-2 moves, which is the most relevant setting for computation. 

Fundamentally this is an experimental paper, though the theoretical underpinnings are in- 
teresting in their own right. Based on an exhaustive census of almost 150 million triangulations, 
including 81800 394 one- vertex triangulations of 1901 distinct 3-manifolds, the answers to the 
questions above appear to be: 

(i) we can relate and simplify one- vertex triangulations using remarkably few Pachner moves, 
and the average number of moves decreases as the number of tetrahedra grows; 

(ii) both procedures require us to add at most two extra tetrahedra, which explains why greedy 
simplification works so well; 

(iii) the number of moves required in the worst case to simplify a 3-sphere triangulation grows 
extremely slowly, to the point where sub-exponential time 3-sphere recognition may indeed 
be possible. 

These observations are extremely surprising, especially in light of Mijatovic's bounds. For 
arbitrary manifolds, observation (ii) does not generalise: in Section [5] we construct larger trian- 
gulations of graph manifolds, beyond the limits of our census, for which three extra tetrahedra 
are required. In the case of the 3-sphere, no such counterexamples are known. If (iii) can be 
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proven in general — yielding a sub-exponential time 3-sphere recognition algorithm — this would 
be a significant breakthrough in computational topology. 

In Section[2]we outline preliminary concepts and introduce Pachner graphs, which are infinite 
graphs whose nodes represent triangulations and whose arcs represent Pachner moves. These 
graphs are the framework on which we build the rest of the paper. We define simplification 
paths through these graphs, as well as the key quantities of length and excess height that we 
seek to measure. 

We follow in Section |3] with two key tools for studying Pachner graphs: an isomorph-free 
census of all closed 3-manifold triangulations with < 9 tetrahedra (which gives us the nodes of 
the graphs), and isomorphism signatures of triangulations that can be computed in polynomial 
time (which allow us to construct the arcs of the graphs). Here we also prove that the census 
grows at a super-exponential rate, despite its strong topological constraints. 

Section d] introduces theoretical techniques for "reducing" paths through Pachner graphs, 
and describes parallel algorithms that bound both the length and excess height of such paths. 
These algorithms are designed to work within the severe time and memory constraints imposed 
by the super-exponential census growth rate. This section also presents the highly unexpected 
experimental results outlined above. 

In Section [5] we study pathological cases, including the census triangulations that are most 
difficult to simplify, as well as the graph manifold constructions mentioned above. We finish in 
Section [6] by exploring the wider implications of our experimental results, in particular for the 
worst-case and generic complexity analysis of topological decision problems. 

All code was written using the topological software package Regina [31 [TU] . 

2 Triangulations and the Pachner graph 

A 3-manifold triangulation of size 7i is a collection of n tetrahedra whose An faces are affinely 
identified (or "glued together") in 2n pairs so that the resulting topological space is a closed 
3-manifold0 We are not interested in the shapes or sizes of tetrahedra (since these do not affect 
the topology), but merely the combinatorics of how the faces are glued together. Throughout 
this paper, all triangulations and 3-manifolds are assumed to be connected. 

We do allow two faces of the same tetrahedron to be identified, and we also note that distinct 
edges or vertices of the same tetrahedron might become identified as a by-product of the face 
gluings. A set of tetrahedron vertices that are identified together is collectively referred to as a 
vertex of the triangulation; we define an edge or face of the triangulation in a similar fashion. 




Figure 1: A 3-manifold triangulation of size n = 2 

^It is sometimes useful to consider bounded triangulations where some faces are left unidentified, or ideal 
triangulations where the overall space only becomes a 3-manifold when we delete the vertices of each tetrahedron. 
Such triangulations do not concern us here. 
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Figure [T] illustrates a 3- manifold triangulation of size n — 2. Here the back two faces of 
the first tetrahedron are identified with a twist, the front faces of the first tetrahedron are 
identified with the front faces of the second using more twists, and the back faces of the second 
tetrahedron are identified together by directly "folding" one onto the other. This is a one-vertex 
triangulation since all eight tetrahedron vertices become identified together. The triangulation 
has three distinct edges, indicated in the diagram by three distinct arrowheads. 

For a given 3-manifold M, a minimal triangulation of is a triangulation of Ai that uses 
the fewest possible tetrahedra. 

Not every pairwise gluing of tetrahedron faces results in a 3-manifold triangulation. Given 
n tetrahedra whose faces are affinely identified in pairs, we obtain a 3-manifold triangulation if 
and only if: (i) every vertex of the triangulation has a small regular neighbourhood bounded by 
a sphere (not some higher-genus surface) , and (ii) no edge of the triangulation is identified with 
itself in reverse. 

The four Pachner moves describe local modifications to a triangulation. These include: 

• the 2-3 move, where we replace two distinct tetrahedra joined along a common face with 
three distinct tetrahedra joined along a common edge; 

• the 1-4 move, where we replace a single tetrahedron with four distinct tetrahedra meeting 
at a common internal vertex; 

• the 3-2 and 4-1 moves, which are inverse to the 2-3 and 1-4 moves. 

These four moves are illustrated in Figure [H Essentially, the 1-4 and 4-1 moves retriangulate 
the interior of a pyramid, and the 2-3 and 3-2 moves retriangulate the interior of a bipyramid. It 
is clear that Pachner moves do not change the topology of the triangulation (i.e., the underlying 
3-manifold remains the same). Another important observation is that the 2-3 and 3-2 moves do 
not change the number of vertices in the triangulation. 



Two triangulations are isomorphic if they are identical up to a relabelling of tetrahedra and a 
reordering of the four vertices of each tetrahedron (that is, isomorphic in the usual combinatorial 
sense). Up to isomorphism, there are finitely many distinct triangulations of any given size. 

Pachner originally showed that any two triangulations of the same closed 3-manifold can be 
made isomorphic by performing a sequence of Pachner moves 134 J Matveev later strengthened 
this result to show that any two one-vertex triangulations of the same closed 3-manifold with at 

■^As Mijatovic notes, Pachner's original result was proven only for true simplicial complexes, but it is easily 
extended to the more flexible deflnition of a triangulation that we use here 1291 . The key step is to remove 
irregularities by performing a second barycentric subdivision using Pachner moves. 




(a) The 2-3 and 3-2 moves 



(b) The 1-4 and 4-1 moves 



Figure 2: The four Pachner moves for a 3-manifold triangulation 
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least two tetrahedra can be made isomorphic through a sequence of 2-3 and/or 3-2 moves 
The two-tetrahedron condition is required because it is impossible to perform a 2-3 or 3-2 move 
upon a one-tetrahedron triangulation (each move requires two or three distinct tetrahedra) . 

In this paper we introduce the Pachner graph, which describes how distinct triangulations 
of a closed 3-manifold can be related via Pachner moves. We define this graph in terms of nodes 
and arcs, to avoid confusion with the vertices and edges that appear in 3-manifold triangulations. 

Definition (Pachner graph). Let AI be any closed 3-manifold. The Pachner graph of M, 
denoted ^(Af ), is an infinite graph constructed as follows. The nodes of ^{AI) correspond to 
isomorphism classes of triangulations of M. Two nodes of ^{M) are joined by an arc if and 
only if there is some Pachner move that converts one class of triangulations into the other. 

The restricted Pachner graph of M, denoted ^i(M), is the subgraph of ^{AI) defined by 
only those nodes corresponding to one-vertex triangulations. The nodes of ^{AI) and £Pi{AI) 
are partitioned into finite levels 1,2,3,..., where each level n contains the nodes corresponding 
to n-tetrahedron triangulations. 




Level 3 
Level 2 

• Level 1 

Figure 3: Levels 1-3 of the restricted Pachner graph of the 3-sphere 

It is clear that the arcs are well-defined (since Pachner moves are preserved under isomor- 
phism), and that arcs do not need to be directed (since each 2-3 or 1-4 move has a corresponding 
inverse 3-2 or 4-1 move). In the full Pachner graph ^{M), each arc runs from some level i to 
a nearby level j ± 1 or i ± 3. In the restricted Pachner graph ^i(M), each arc must describe a 
2-3 or 3-2 move, and must run from some level i to an adjacent level i±l. Figure |3] shows the 
first few levels of the restricted Pachner graph of the 3-sphere. 

We can now reformulate the results of Pachner and Matveev as follows: 

Theorem 2.1 (Pachner, Matveev). The Pachner graph of any closed 3-manifold is connected. 
If we delete level 1, the restricted Pachner graph of any closed 3-manifold is also connected. 

To simplify a triangulation we essentially follow a path through ^(M) or ^i{AI) from a 
higher level to a lower level, motivating the following definitions: 

Definition (Simplification path). A simplification path is a directed path through either ^^(M) 
or 3^i{AI) from a node at some level i to a node at some lower level < i. 

Definition (Length and excess height). Let p be any path through ^(M) or 3^i{AI) from 
level i to level j. The length of p is the number of arcs it contains. The excess height of p is the 
smallest /i > for which the entire path stays in or below level (max{i, j} -I- h). 

Figure ID illustrates a path of length 13 and excess height 2. For simplification paths, the 
length and excess height measure how difficult it is to simplify a triangulation: the length 
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Level 9 




T 



Level 8 



excess height = 2 



Level 7 = i 



Level 6 = j 
Level 5 



Figure 4: Measuring length and excess height 



measures the number of Pachner moves, and the excess height measures the number of extra 
tetrahedra required. 

To simplify a 3-sphere triangulation, the only known bounds on length and excess height are 
due to Mijatovic [29!: 

Theorem 2.2 (Mijatovic). Any triangulation of the 3-sphere can be converted into a two- 
tetrahedron triangulation using less than 6 • l0'^n'^2'^-^° " Pachner moves. 

Corollary 2.3. In the Pachner graph ^{S^), from any node at level n > 2 there is a simplifi- 
cation path of length less than 6 • 10^n^2^'^'^ " and excess height less than 3 • 10^n^2^'^'^ " . 

Mijatovic also proves bounds for other classes of 3-manifolds, including Seifert fibred spaces 
[50] . fibre- free Haken manifolds [32] j and knot complements [31- These all involve tovifers of 
exponentials (where in some cases the height of the tower grows with n), and the resulting 
bounds are far greater than the 3-sphere bounds cited above. 

In the restricted Pachner graph, where we only consider 2-3 and 3-2 moves, no explicit 
bounds on length or excess height are known for any 3-manifolds at all. 

3 Key tools 

Experimental studies of Pachner graphs are difficult: the graphs themselves are infinite, and 
even the finite level sets grow super-exponentially in size (as we show in Theorem 13. ip . By 
working with isomorphism classes of triangulations, we keep the level sets considerably smaller 
than if we had used labelled triangulations instead. However, the trade-off is that both the 
nodes and the arcs of each graph become more difficult to construct. 

In this section we outline two key algorithmic tools for studying Pachner graphs: a census 
of triangulations (which enumerates the nodes at each level), and polynomial-time computable 
isomorphism signatures (which allow us to construct the arcs). 

3.1 A census of triangulations 

To enumerate the nodes of Pachner graphs, we build an exhaustive census of all 3-manifold 
triangulations of size n <9, with each triangulation included precisely once up to isomorphism. 
The total number of triangulations in this census is 149 676 922. Following the focus of this 
paper, we extract from these the one-vertex triangulations of prime orientable 3-manifolds and 
the 3-sphere. The resulting 81 800 394 triangulations represent 1 901 distinct 3-manifolds, and 
these triangulations form the basis of our experiments. 
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107 185 448 
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4472 


50 778 389 










81 800 394 of interest 



Table 1: Breakdown of different types of triangulations in the census 



A full breakdown of triangulations in the census appears in Table [T] For prime orientable 
3-manifolds, we further divide these triangulations into non-minimal (where we study simpli- 
fication paths) versus minimal (where we study how difficult it is to join distinct minimal 
triangulations of the same 3-manifold). 

The algorithms behind this census are sophisticated: 

• Generating triangulations: This is a combinatorial enumeration problem with severe topo- 
logical constraints. If we simply enumerate all pairwise identifications of tetrahedron faces 
up to isomorphism, there are at least 

[(4n - 1) X (4n - 3) X ■ ■ ■ X 3 X 1] ■ 6^" ^ ^ 35 x IQ^^ 
n\ ■ 24" 

possibilities for n = 9. However, the topological constraint that the triangulation must 
represent a 3-manifold cuts this number down to just ~ 1.39 x 10^, as seen in Table [TJ 

A key challenge therefore is to enforce this topological constraint as the census runs, and 
thus prune vast branches of the combinatorial search tree. Techniques for this include 
modified union-find and skip list algorithms for tracking partially-constructed edge and 
vertex links ^ [S] , and the analysis of 4-valent face pairing graphs [H d] . Some authors 
describe other techniques specific to minimal triangulations [231 HSl US] ) but these are too 
specialised for the larger body of data that we require here. 

For the largest case n = 9, the full enumeration of triangulations required ^ 85 days of 
CPU time as measured on a single 1.7 GHz IBM Power5 processor. In reality this was 
reduced to 2-3 days of wall time using 32 CPUs in parallel. 

The paper [8j expands this census to n = 10 (with 2 046 869 999 triangulations), but we 
do not use this data here because the resulting Pachner graphs are too large to process. 
See Section |4] for further discussion on time and memory constraints. 



Identifying 3-spheres and orientable prime 3-manifolds: Both 3-sphere recognition and 
prime decomposition are theoretically difficult problems. Although the best known algo- 
rithms have exponential running times, recent advances have made them extremely fast 
for problems of our size. 
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We employ algorithms based on the 0-efficiency techniques of Jaco and Rubinstein [T5] . 
coupled with highly optimised algorithms for normal surface enumeration [B] and 3-sphere 
recognition [7]. The total running time over all ~ 150 million triangulations was just 
7.7 hours, which is negligible in comparison to the census enumeration. 

• Identifying minimal and non-minimal triangulations: For this we call upon a separate, 
specialised census of minimal 3-manifold triangulations. Minimal triangulations are ex- 
tremely rare (as seen Table [1]), which means that there are more opportunities for pruning 
the combinatorial search tree, and so specialised censuses of minimal triangulations are 
significantly faster to build. 

Censuses with at least one minimal triangulation per prime orientable 3-manifold have 
been compiled by Matveev and others for n < 12 tetrahedra [35], and censuses of all 
minimal triangulations of such manifolds have been compiled by the author for n < 11 
tetrahedra J^. Since this study requires all minimal triangulations, we use [5J as our 
source. 

It was noted in the first point above that 3-manifold triangulations are extremely rare 
amongst all pairwise identifications of tetrahedron faces. Dunfield and Thurston [T2] justify 
this theoretically, proving that as n cx), the probability that a random identification of tetra- 
hedron faces yields a 3-manifold triangulation tends to zero. Despite this rarity, we can prove 
that our census of 3-manifold triangulations grows at a super-exponential rate: 

Theorem 3.1. The number of distinct isomorphism classes of 3-manifold triangulations of size 
n grows at an asymptotic rate o/exp(0(nlog7i)). 

Proof. An upper bound of exp(0(n log n)) is easy to obtain. If wc count all possible pairwise 
gluings of tetrahedron faces, without regard for isomorphism classes or other constraints (such 
as the need for the triangulation to represent a closed 3-manifold), we obtain an upper bound 
of 

[{An - 1) X (4n - 3) X • • • X 3 X 1] • 6^" < (4^)^" • 6^" £ exp(0(nlog7i)). 

Proving a lower bound of exp(f2(7ilog?i)) is more difficult — the main complication, as noted 
above, is that most pairwise identifications of tetrahedron faces do not yield a 3-manifold at 
all. We work around this by first counting 2-manifold triangulations (which are easier to ob- 
tain), and then giving a construction that "fattens" these into 3-manifold triangulations without 
introducing any unwanted isomorphisms. 

To create a 2-manifold triangulation of size 2m (the size must always be even), we identify 
the 6m edges of 2m distinct triangles in pairs. Any such identification will always yield a closed 
2-manifold; that is, nothing can "go wrong" , in contrast to the three-dimensional caseH 

There is, however, the issue of connectedness to deal with (recall from the beginning of 
Section [2] that all triangulations in this paper are assumed to be connected) . To ensure that a 
labelled 2-manifold triangulation is connected, we insist that for each fc = 2, 3, . . . , 2m, the first 
edge of the triangle labelled k is identified with some edge from one of the triangles labelled 
1, 2, . . . , A; — 1. Of course many connected labelled 2-manifold triangulations do not have this 
property, but since we are proving a lower bound this does not matter. 

We can now place a lower bound on the number of labelled 2-manifold triangulations. First 
we choose which edges to pair with the first edges from triangles 2,3,..., 2m; from the property 

•^Recall that things "go wrong" in three dimensions if a small neighbourhood of a vertex is not surrounded 
by a sphere, or if an edge is identified with itself in reverse. 
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above we have 3 x 4 x . . . x 2m x (2m + 1) = ^ (2m + 1) ! choices. We then pair off the remaining 
2m + 2 edges, with (2m + 1) x (2m — l)x...x3xl = (2m + l)!/2™m! possibihties overall. 
Finally we note that each of the 3m pairs of edges can be identified using one of two possible 
orientations. The total number of labelled 2-manifold triangulations is therefore at least 

(2m +1)! (2m + 1)! _ (2m + 1)! ■ (2m + 1)! • 22" 

2 2™m! ' ^ 2 • m! ' 

Each isomorphism class contains at most (2m)! • 6^™ labelled triangulations, and so the 
number of distinct isomorphism classes of 2-manifold triangulations is bounded below by 

(2m + 1)! ■ (2m + 1)! • 2^" _ (2m + 1) • (2m + 1)! 
2 • m! • (2m)! • G^™ ~ 2 • m! • 3^™ 

> (2m + 1) X 2m X • • • X (m + 2) X (m + 1) X (g)"' 
>(m + l)"+i.(i)" 
G exp(r2(mlogm)). 

We fatten each 2-manifold triangulation into a 3-manifold triangulation as follows. Let F 
denote the closed 2-manifold described by the original triangulation. 

1. Replace each triangle with a prism and glue the vertical faces of adjacent prisms together, 



as illustrated in Figure 5(a) This represents a bounded 3-manifold, which is the product 
space F X I. 



2. Cap each prism at both ends with a triangular pillow, as illustrated in Figure 5(b) The 
two faces of each pillow are glued to the top and bottom of the corresponding prism, 
effectively converting each prism into a solid torus. This produces the closed 3-manifold 



F X S^, and the complete construction is illustrated in Figure 5(c) 



3. Triangulate each pillow using two tetrahedra, which are joined along three internal faces 
surrounding an internal vertex. Triangulate each prism using 14 tetrahedra, which again 
all meet at an internal vertex. Both triangulations are illustrated in Figure [5 (d)[ 

If the original 2-manifold triangulation uses 2m triangles, the resulting 3-manifold triangu- 
lation uses n = 32m tetrahedra. Moreover, if two 3-manifold triangulations obtained using this 
construction are isomorphic, the original 2-manifold triangulations must also be isomorphic. 
The reason for this is as follows: 

• Any isomorphism between two such 3-manifold triangulations must map triangular pillows 
to triangular pillows. This is because the internal vertex of each triangular pillow meets 
only two tetrahedra, and no other vertices under our construction have this property. 

• By "flattening" the triangular pillows into 2-dimensional triangles, we thereby obtain an 
isomorphism between the underlying 2-manifold triangulations. 

It follows that, for n = 32m, we obtain a family of exp(i7(mlogm)) = exp(i7(n log n)) 
pairwise non-isomorphic 3-manifold triangulations. 

This result is easily extended to n ^ mod 32. Let Vn denote the number of distinct 
isomorphism classes of 3-manifold triangulations of size n. 
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(c) The complete construction (d) Triangulating prisms and pillows 



Figure 5: Fattening a 2-nianifold triangulation into a 3- manifold triangulation 

• Each triangulation of size n has at least n — 1 distinct 2-3 moves available (since any face 
joining two distinct tetrahedra defines a 2-3 move, and there are at least n — 1 such faces). 

• On the other hand, each triangulation of size n+1 has at most 6(n-|- 1) distinct 3-2 moves 
available (since each 3-2 move is defined by an edge that meets three distinct tetrahedra, 
and the triangulation has at most 6(n -t- 1) edges in total). 

It follows that Vn+i > Vn ■ > Vn/18 for any n > 1. This gives V32,n+k > V32m/18^^ 

for sufficiently large m and all < fc < 32, and so we obtain Vn G exp(f2(nlogn)) with no 
restrictions on n. □ 

Remark. Of course, we expect that Vn+i ^ Vn (and indeed we see this in the census). The 
bounds that we use to show Vi+i > Vn/18 in the proof above are very loose, but they are 
sufiicient for the asymptotic result that we seek. 

3.2 Isomorphism signatures 

To construct the arcs of a Pachner graph, we begin at a node — that is, a 3-manifold triangulation 
T — and perform Pachner moves. Our main difficulty is determining the endpoints of the arcs: 
each Pachner move results in a new triangulation 7"', and we must determine which node of the 
graph represents T' ■ 

A naive approach might be to search through all nodes at the appropriate level of the Pachner 
graph and test each triangulation for isomorphism with T' . However, this is infeasible: although 
isomorphism testing is fast (as we prove in Corollary 13. 5p . the sheer number of nodes at level n 
of the graph is too large (as shown by Theorem 13. 

What we need is a property of the triangulation T' that is easy to compute, and that uniquely 
defines the isomorphism class of T' . This property can be used as the key in a data structure 
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with fast insertion and fast lookup (such as a hash table or a red-black tree), and by computing 
this property we can quickly jump to the relevant node of the graph. 

Here we define such a property, which we call the isomorphism signature of a triangulation. 
In Theorem 13 . 31 we show that isomorphism signatures uniquely define isomorphism classes, and 
in Theorem 13.41 we show that they are small to store and fast to compute. 

A labelling of a triangulation of size n involves: (i) numbering its tetrahedra from to n — 1 
inclusive, and (ii) numbering the four vertices of each tetrahedron from to 3 inclusive0 We 
also label the four faces of each tetrahedron from to 3 inclusive so that face i is opposite vertex 
i. A key ingredient of isomorphism signatures is canonical labellings, which we define as follows. 

Definition (Canonical labelling). Given a labelling of a triangulation of size n, let Atj denote 
the tetrahedron glued to face / of tetrahedron t (so that Atj £ {0, ...,n — 1} for all t — 
0, . . . , n — 1 and / = 0, . . . , 3). The labelling is canonical if, when we write out the sequence 
^0,0, ^0,1, ^0,2, ^0,3, ^1,0, ■ • ■ , the following properties hold: 

(i) For each 1 < i < j, tetrahedron i first appears before tetrahedron j first appears. 

(ii) For each i > 1, suppose tetrahedron i first appears as the entry Atj = i. Then the 
corresponding gluing uses the identity map: face / of tetrahedron t is glued to face / of 
tetrahedron i so that vertex v of tetrahedron t maps to vertex v of tetrahedron i for each 

As an example, consider the triangulation of size n — 3 described by Table H This table 
lists the precise gluings of tetrahedron faces. For instance, the second cell in the bottom row 
indicates that face 1 of tetrahedron 2 is glued to tetrahedron 1, in such a way that vertices 0, 2, 3 
of tetrahedron 2 map to vertices 3,1,2 of tetrahedron 1 respectively. This same gluing can be 
seen from the other direction by examining the first cell in the middle row. 





Face 
Vertices 123 


Face 1 
Vertices 023 


Face 2 
Vertices 013 


Face 3 
Vertices 012 


Tet. 
Tet. 1 
Tet. 2 


Tet. 


120 


Tet. 1 
Tet. 


023 
023 


Tet. 2 
Tet. 1 
Tet. 


013 
012 
013 


Tet. 
Tet. 1 
Tet. 2 


312 
013 
123 


Tet. 2 


230 


Tet. 2 


012 


Tet. 1 


312 



Table 2: The tetrahedron face gluings for an example 3-tetrahedron triangulation 

It is simple to see that the labelling for this triangulation is canonical. The sequence 
^o,Oj • ■ • , is 0,1,2,0,2,0,1,1,2,1,0,2 (reading tetrahedron numbers from left to right 

and then top to bottom in the table), and tetrahedron 1 first appears before tetrahedron 2 
as required. Looking closer, the first appearance of tetrahedron 1 is in the second cell of the 
top row where vertices 0, 2, 3 map to 0, 2, 3, and the first appearance of tetrahedron 2 is in the 
subsequent cell where vertices 0,1,3 map to 0, 1, 3. In both cases the gluings use the identity 
map. 

Lemma 3.2. For any triangulation T of size n, there are precisely 24n canonical labellings of 
T, and these can be enumerated in O(n^logn) time. 

*We start numbering from instead of 1 for consistency with the software implementation of isomorphism 
signatures, as detailed in the appendix. 
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Proof. For n = 1 the result is trivial, since all 24 — 4! possible labellings are canonical. For 
n > 1 we observe that, if we choose (i) any one of the n tetrahedra to label as tetrahedron 0, 
and (ii) any one of the 24 possible labellings of its four vertices, then there is one and only one 
way to extend these choices to a canonical labelling of T. 

To see this, we can walk through the list of faces Fq^q, Fo,i^Fo,2, -Po,3i -P'l.Oi • ■ • i Fn-i,3, where 
Ft i represents face i of tetrahedron t. The first face amongst i^o,Oj • ■ • i ^o,3 that is joined to an 
unlabelled tetrahedron must in fact be joined to tetrahedron 1 using the identity map. This 
allows us to deduce tetrahedron 1 as well as the labels of its four vertices. 

We inductively extend the labelling in this manner: once we have labelled tetrahedra 0, . . . , fc 
and their corresponding vertices, the first face amongst Fq^, . . . , Fk^3 that is joined to an un- 
labelled tetrahedron must give us tetrahedron fc + 1 along with the labels for its four vertices 
(again using the identity map). The resulting labelling is canonical, and all of the labels can 
be deduced in 0{n\ogn) time using a single pass through the list Fo_o, . • . , i^„-i,3. The logn 
factor is for manipulating tetrahedron labels, each of which requires O(logn) bits. 

It follows that there are precisely 24n canonical labellings of T, and that these can be 
enumerated in 0(n^ log n) time using 24n iterations of the procedure described above. □ 

Definition (Isomorphism signature). For any triangulation T of size n, enumerate all 24n 
canonical labellings of T, and for each canonical labelling encode the full set of face gluings as 
a sequence of bits. We define the isomorphism signature to be the lexicographically smallest of 
these 24n bit sequences, and we denote this by (j{T). 

For the theoretical results in this section, it suffices to treat each sequence as a naive bitwise 
encoding of a full table of face gluings (such as Table [5]). In practical software settings however, 
we use a more compact alphanumeric representation with less redundancy, which we specify in 
full detail in the appendix. 

Theorem 3.3. Given two 3-manifold triangulations T and T' , we have cr(T) — c(7~') if and 
only if T and T' are isomorphic. 

Proof. It is clear that (j{T) — cr{T') implies that T and T' are isomorphic, since both signatures 
encode the same gluing data. Conversely, if T and T' are isomorphic then their 24n canonical 
labellings are the same (though they might be enumerated in a different order). In particular, 
the lexicographically smallest canonical labellings will be identical; that is, (t(T) = ct(T')- D 

Theorem 3.4. Given a 3-manifold triangulation T of size n, the isomorphism signature (t(7~) 
has 0(n log n) size and can be generated in 0{n^\ogn) time. 

Proof. To encode a full set of face gluings, at worst we require a table of gluing data such as 
Table [21 with 4n cells each containing four integers. Because some of these integers require 
O(logn) bits (the tetrahedron labels), it follows that the total size of cr(T) is O(rilogri). 

The algorithm to generate a{T) is spelled out explicitly in its definition. The 24n canonical 
labellings of T can be enumerated in 0{n^ logn) time (Lemma 13. 2|) . Because a full set of face 
gluings has size O(rilogn), we can encode the 24n bit sequences and select the lexicographically 
smallest in 0{n^ logn) time, giving a time complexity of 0{n^ \ogn) overall. □ 

This space complexity of 0{n log n) is the best we can hope for, since Theorem 13 . 1 1 shows that 
the number of distinct isomorphism signatures for size n triangulations grows like exp(0(n log n)) . 

It follows from Theorems 13.31 and 13.41 that isomorphism signatures are ideal tools for con- 
structing arcs in the Pachner graph, as explained at the beginning of this section. Moreover, 
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the relevant definitions and results are easily extended to bounded and ideal triangulations (as 
described in the appendix). We finish with a simple but important consequence of our results: 

Corollary 3.5. Given two 3-manifold triangulations T and T' each of size n, we can test 
whether T and T' are isomorphic in 0{n^ logn) time. 

4 Analysing Pachner graphs 

As discussed in the introduction, our focus is on one-vertex triangulations of closed prime 
orientable 3-manifolds and of the 3-sphere. We therefore direct our attention to the restricted 
Pachner graph 3^i{AA) of each such 3-manifold AA. 

Definition (Base level). For any closed 3-manifold A^, the base level /3i(A^) is the lowest 
non-empty level of the restricted Pachner graph ^i(A^), excluding the isolated level 1. In 
other words, j3i{M.) is the smallest n > 2 for which there exists an n-tetrahedron, one- vertex 
triangulation of M.. 

Lemma 4.1. For any closed prime orientable 3-manifold Ai except for the lens spaces -/j(4, 1) 
and L(5,2), the base level j3i{A4) is simply the size of a minimal triangulation of M. If M = 
L(4, 1) or L(5,2), then l3i{M) — 3. For the 3-sphere (which is not prime), f3i{S'^) = 2. 

Proof. From Table [H there are only three one- vertex triangulations of closed 3-manifolds; these 
represent , L(4, 1) and L(5, 2). The next smallest one- vertex triangulations of these manifolds 
in the census have sizes fii{S^) = 2, ^\(L{^, 1)) = 3 and /3i(L(5, 2)) = 3. 

For any other closed prime orientable 3-manifold A^, either M. = MP'^, M. = L(3,l), or 
every minimal triangulation of has one vertex [25]. In the latter case, the result follows 
immediately from the definition of /3i(A^). For both MP^ and L(3, 1) at least one minimal 
triangulation has one vertex, and so the result holds for these cases also. □ 

In this section, we analyse all 81 800 394 triangulations of interest of size n < 9 in our census 
to obtain computer proofs of the following results: 

Theorem 4.2 (Excess heights of simplification paths). Let be a closed prime orientable 
3-manifold or the 3-sphere. From any node of ^i{M) at level n where /3i{M) < n < 9, there 
is a simplification path of excess height < 2. 

Theorem 4.3 (Lengths of simplification paths). Let M be a closed prime orientable 3-manifold. 
From any node of at level n where ^i(Al) < n < 9, there is a simplification path of 

length < 17. For the 3-sphere, the bound is stronger: from any node of 3^i{S^) at level n where 
PiiS^) <n<9, there is some simplification path of length < 9. 

Theorem 4.4 (Joining minimal triangulations). Let A4 be a closed prime orientable 3-manifold 
that is not the lens space L(3, 1), and for which /3i(A4) < 9. Then any two nodes at level /3i(A4) 
of ^i(A^) are joined by a path of length < 18 and excess height < 2. 

For the special case L(3, 1), there are precisely two nodes at level /?i(L(3, 1)) = 2, and these 
are joined by a path of length 6 and excess height 3. 

These results are astonishing, given Mijatovic's super-exponential bounds for the 3-sphere 
P5] and tower-of-exponential bounds for other manifolds [5D1 [3TJ [32] . Theorems 14.21 and 14.41 
state that we can simplify triangulations or convert between minimal triangulations using at 
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most two extra tetrahedra, and Theorems 14.31 and 14.41 state that we require very few Pachner 
moves. The 3-sphere resuhs in particular have important algorithmic implications, which we 
discuss further in Section [51 

The algorithms behind these computer proofs are specialised, and use novel techniques to 
control the enormous time and space requirements. We return to these algorithmic issues shortly. 
In the meantime, some further comments on Theorems I4.2| - H^ 

• We restrict our results to levels n < 9 because, even though we have census data for n = 10 
(see [5]), the sheer size of the census makes it infeasible to compute the necessary properties 
of the Pachner graphs. As a result. Theorems 14.21 and 14.31 cover 747 distinct 3-manifolds 
with (3i{M) < 9, and Theorem li^ covers 1900 distinct 3-manifolds with (3i{M) < 9. 

• Theorem |43] bounds the worst-case length of the shortest simplification path. If we bound 
the average length instead, the results are much smaller still: for the case n = 9, this 
length is less than 1.89 when averaged over all triangulations of S^, and less than 1.91 
when averaged over all triangulations of closed prime orientable 3-manifolds. We present 
these average-case results in detail in Section 14.31 

• For arbitrary closed prime orientable 3-manifolds, the height results do not generalise. In 
Section [SJ we construct (i) a triangulation of a graph manifold of size n = 10 for which 
every simplification path has excess height > 3, and (ii) two minimal triangulations of a 
graph manifold with size n = 10 where every path joining them has excess height > 3. 
For the case of the 3-sphere, no such counterexamples are known. 

For the remainder of this section, we describe the algorithms behind Theorems I4.2H4.41 and 
present the experimental results in detail. Our algorithms are constrained by the following 
factors: 

• Their time and space complexities must be close to linear in the number of nodes that 
they examine, due to the sheer size of the census. 

• They cannot loop through all nodes in ^i{Ai), since the graph is infinite. They cannot 
even loop through all nodes at any level n > 11, since there are too many to enumerate. 

• They cannot follow arbitrary breadth- first or depth- first searches through ^i{A4), since 
the graph is infinite and can branch heavily in the upward directionlf] 

Because of these limiting factors, we cannot run through the census and directly measure 
the shortest length or smallest excess height of any simplification path from each node. Instead 
we develop fast, localised algorithms that allow us to bound these quantities from above. To 
our delight, these upper bounds are extremely effective in practice. 

A key optimisation in many of our algorithms comes from Theorem 14.51 which allows us to 
"reduce" sequences of Pachner moves to use fewer intermediate tetrahedra. We state and prove 
this theorem in Section [4.1l Following this, we present the individual algorithms: 

• In Section 14.21 we give a fast algorithm based on union-find that bounds the heights of 
simplification paths. We run this algorithm over the census of 3-sphere triangulations, 
giving a computer proof for the 3-sphere case of Theorem 14.21 

^In general, a node at level n can have up to 2n distinct neighbours at level (ra -|- 1). 
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• In Section 14.31 we describe a multiple-source breadth-first search algorithm that bounds 
lengths of simplification paths. By running this over the census of triangulations of the 
3-sphere and closed prime orientable manifolds, we prove Theorem 14.31 as well as the 
outstanding closed prime orientable case of Theorem 14.21 

• In Section IT^ we present algorithms that incorporate both breadth-first search and union- 
find techniques to study paths between nodes representing minimal triangulations. By 
running these algorithms over the census of minimal triangulations of closed prime ori- 
entable manifolds, we obtain a computer proof of Theorem 14.41 

• We finish in Section 14.51 with a discussion of how these algorithms can be parallelised 
effectively, along with explicit measurements of running time and memory use. 

4.1 Reducing sequences of moves 

The key idea behind reducing sequences of Pachner moves is that, in many cases, we can 
interchange pairs of consecutive 2-3 and 3-2 moves to become consecutive 3-2 and 2-3 moves 
instead, without changing the final triangulation. In the Pachner graph, this replaces a path 
segment from levels n — >■ (n -I- 1) — > n with a path segment from levels n — > (n — 1) — > n: the 
length and endpoints of the overall path remain the same, and the excess height will either stay 
the same or decrease (depending on the heights of other segments of the path) . 

This type of interchange is not always possible. The purpose of Theorem 14.51 is to identify 
the "bad cases" where consecutive 2-3 and 3-2 moves cannot be interchanged, so that we can 
explicitly test for them in algorithms. In summary, we find that every bad 2-3 / 3-2 pair can be 
described by one of three "composite moves" . These composite moves are local modifications to 
a triangulation, much like the Pachner moves (though a little more complex), and are defined 
as follows. 

Definition (Octahedron flip). An octahedron flip, also known as a 4-4 move, involves four 
distinct tetrahedra surrounding a common edge of degree four. The move essentially rotates the 
configuration, replacing it with four new tetrahedra surrounding a new edge of degree four that 
points in a different direction, as illustrated in Figure [6] 



Figure 6: An octahedron flip 



Definition (Pillow flip). 



A quadrilateral pillow is formed from two distinct tetrahedra joined 

A pillow flip is a move involving three 



along two adjacent faces, as shown in Figure 7(a) 
tetrahedra: a quadrilateral pillow plus a third tetrahedron attached to one of its outer faces. 
The move essentially reflects the conflguration, replacing it with a quadrilateral pillow plus a 
new tetrahedron attached to the opposite face instead. 
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(a) A quadrilateral pillow (b) Performing a pillow flip 

Figure 7: A pillow flip 



This move is illustrated in Figure [7(b)[ in which the pillows are shaded. In the left diagram 
the pillow is at the front, and the third tetrahedron is attached to its lower rear face. In the 
right diagram the pillow moves to the rear, and the third tetrahedron is attached to its lower 
front face instead. 

Definition (Prism flip). There are two types of prism flip, which we call types A and B. 
Both begin with three distinct tetrahedra joined to form a triangular prism, as illustrated in 
Figure 8(a) For the type A flip we also require that two rectangular faces of the prism are 
folded together, as shown in Figure [8(b)] for type B we require that these same two rectangular 
faces be folded together with a 180° twist instead. In both cases, the prism flip involves rotating 
the entire configuration so that the two triangular ends of the prism are interchanged, as shown 



in Figure 8(c) 



(a) A 3-tetrahedron triangular prism (b) Folding two rectangular faces together 




(c) Performing the move 

Figure 8: A prism flip of type A 



Remark. It is easy to search for all possible octahedron, pillow or prism flips on a given 
triangulation, since each occurs around an edge of degree four, two or five respectively. We can 
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simply search for edges of the correct degree, test for the necessary configuration of tetrahedra, 
and perform the flip if possible. 

Note that there could be two different octahedron flips around the same degree four edge 
(corresponding to the two possible directions for the new edge that replaces it), and there could 
be four different pillow flips around the same degree two edge (corresponding to the four faces 
of the pillow to which the third tetrahedron might be attached). Around a degree five edge, 
there can only ever be one prism flip (if there is any at all). 

Theorem 4.5. Let T and T' he 3-manifold triangulations, where T' can be obtained from T 
by performing a 2-3 move followed by a 3-2 move. Then one of the following cases holds: 

• T' is isomorphic to T; 

• T' can be obtained from T by performing a 3-2 move followed by a 2-3 move instead (i.e., 
using one fewer tetrahedron at the intermediate stage instead of one more ); 

• 7"' can be obtained from T by performing either an octahedron flip, a pillow flip, or a 
prism flip. 

Proof. Suppose that we obtain T' from 7" by: 

(i) performing a 2-3 move on two adjacent tetrahedra Ai, A2 of T, giving the intermediate 
triangulation I; 

(ii) then performing a 3-2 move around the degree three edge e of I to give T' . 

If e is not an edge of the original triangulation 7", then it must be created by the 2-3 
move. Therefore the subsequent 3-2 move is the inverse of the original 2-3 move, and the flnal 
triangulation I'' is isomorphic to T- 

Suppose then that e does belong to the original triangulation T, but that e is not an edge 
of either Ai or A2. This implies that the 2-3 move and the 3-2 move occur in "disjoint" regions 
of the triangulation: none of the tetrahedra involved in the 3-2 move are also involved in the 
2-3 move. Therefore the 2-3 and 3-2 moves can be applied in either order, and we can obtain 
T' from T by performing the 3-2 move followed by the 2-3 move instead. 

Finally, suppose that e is an edge of Ai and/or A2 in T. Note that e might not have degree 
three in T , and it might appear as multiple edges of Ai and/or A2. However, in the intermediate 
triangulation X (after the initial 2-3 move), it must have degree three and it must belong to 
three distinct tetrahedra. 

Up to isomorphism, there are 13 possible ways that e can appear in Ai and/or A2 subject 
to these constraints. These fall into five basic patterns, as illustrated in Figure [3] Cases 1, 2 
and 4 require additional tetrahedra in order to meet the degree three requirement (these extra 
tetrahedra are also pictured). Cases 2, 3 and 5 each have several different subcases, according 
to the specific orientations of e in each tetrahedron and the different ways in which tetrahedron 
faces can be glued together. 

Cases 1 and 4 are simple to analyse. Figure [TU] shows the corresponding configurations in 
T with the extra tetrahedra glued in: these are the initial configurations for an octahedron flip 
and a pillow flip respectively. By following the details of the 2-3 move on Ai and A2 followed 
by the 3-2 move around e, we indeed find that these moves perform a single octahedron flip in 
case 1, and a single pillow fiip in case 4. 

For case 2, there are four possible configurations up to isomorphism. These are shown in 
Figure[TTJ again with the extra tetrahedron glued in. These four subcases correspond to different 
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Case 4 Case 5 

(2 subcases) 



Figure 9: Possible locations for the edge e in T 




Case 1 Case 4 



Figure 10: The full configurations in T for cases 1 and 4 

possible orientations of e in each tetrahedron (indicated in the diagram by the bold arrowheads) 
and different possible choices for which faces are glued together (indicated by the white arrows 
and the shading on the faces). 

Subcase 2a is precisely the initial configuration for a prism flip of type A (the two shaded 
triangles on the left form one rectangular face of the prism, and the two shaded triangles on the 
right form another). Likewise, subcase 2b is the initial configuration for a prism flip of type B 
(the different gluings describe a 180° twist before folding these rectangles together). In each 
case, the 2-3 and 3-2 move together carry out the corresponding prism flip. 

Subcases 2c and 2d both produce non-orientable triangulations. More importantly, each 
contains a vertex V with a non-orientable link — in other words, any small regular neighbourhood 
of V is non-orientable. For subcase 2c, this vertex is at the top of the diagram, and for subcase 2d 
it is the vertex at the centre of the diagram. No such vertex can appear in a 3-manifold 
triangulation, and so subcases 2c and 2d can never occur. 

In case 3, all eight faces of Ai and A2 are glued together in pairs. Since the triangulation T 
is connected, this means that T contains only these two tetrahedra: if n > 2 then this case can 
never occur at all. Even if n = 2, the triangulation is not of interest to us, since T will contain 
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Case 2a Case 2b Case 2c Case 2d 

Figure 11: The four possible configurations for case 2 



two distinct vertices and so will not feature in any restricted Pachner graph. 

Nevertheless, we summarise case 3 with n = 2 for completeness. Here there are five subcases 
up to isomorphism, again depending on the possible orientations of e and choices for how the the 
tetrahedron faces are glued together. Three of these subcases give non-orientable vertex links 
as before, and so cannot occur. Of the final two subcases, one gives T as the unique two-vertex, 
two-tetrahedron triangulation of , and the other gives the unique two- vertex, two-tetrahedron 
triangulation of the lens space L(3, 1). In both subcases, the 2-3 and 3-2 move together produce 
a triangulation that is isomorphic to the original^ 

This leaves case 5. Here there are two possible configurations in T up to isomorphism, as 
illustrated in Figure[T2](in both diagrams, the top left face at the front is glued to the upper face 
at the rear, and the top right face at the front is glued to the lower face at the rear). Subcase 5a 
gives an orientable triangulation (to be precise, a (1, 3, 4) layered solid torus fT9^), and following 
through the 2-3 and 3-2 moves shows once again that the resulting triangulation is isomorphic 
to the original. Subcase 5b gives a non-orientable triangulation in which the top left edge is 
identified with itself in reverse, and so this subcase can never occur. □ 




Case 5a Case 5b 



Figure 12: The two possible configurations for case 5 



4.2 Bounding excess heights 

Our first algorithm computes bounds Hn{A4) for a given 3-manifold M. and a given level n > 
/3i(A^) so that, from every node at level n of the graph there is a simplification path 

of excess height < Hn{M.). By running this algorithm, we explicitly compute these bounds for 
the case Ai = S^, for each n in the range 2 < n < 9. 

^This can be seen because the 2-3 and 3-2 moves together produce another two-vertex, two-tetrahedron 
triangulation of the same manifold, and the census shows that and L{3, 1) each have only one such triangulation 
up to isomorphism. Of course one could simply follow through the moves by hand instead. 
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The general idea is to begin with the set of aU nodes at level n, and repeatedly expand 
into higher levels using 2-3 moves only until all of the original nodes are connected together. 
If the bound Hn{M) is tight (as we find in practice for the 3-sphere), this method of "upward 
expansion" saves significant time and memory by enumerating only a fraction of the higher 
levels £ > n (each of which contains far more nodes than the initial level n). 

Algorithm 4.6 (Computing This algorithm runs by progressively building a finite subgraph 

G C ^i{M.). At all times we keep track of the number of distinct components of G (denoted by 
c) and the maximum level of any node in G (denoted by i). 

1. Initialise G to all of level n of S^i{M.). This means that G has no arcs, the number of 
components c is just the number of nodes at level n, and the maximum level is £ ~ n. 

2. While c > 1, expand the graph as follows: 

(a) Construct all arcs from nodes in G at level i to (possibly new) nodes in ^i{M) at 
level i -\- 1. Insert these arcs and their endpoints into G. 

(b) Update the number of components c, and increment £ by one. 

3. Once we have c = 1, output the final bound Hn{Ai) — £ ~ n and terminate. 

In step 2(a) we construct arcs by performing 2-3 moves, and in step 2(b) we use union-find 
to update the number of components in small time complexity. 

It is clear that Algorithm I4.6l is correct for any n > f3i{A4): once we have c = 1 the subgraph 
G is connected, which means it contains a path from any node at level n to any other node at 
level n. Because n > /3i{Ai), Theorem 12.11 shows that at least one such node allows a 3-2 move, 
and so any node at level n has a simplification path of excess height < £ — n. 

However, it is not clear that Algorithm 14.61 terminates: it might be that every simplification 
path from some node at level n passes through nodes that we never construct at higher levels 
£ > n. Happily it does terminate for the 3-sphere for all 2 < rt < 9, giving an output of 
Hn{S^) = 2 each time, and thereby proving the 3-sphere case for Theorem 14.21 Table [3] shows 
how the number of components c changes as the algorithm runs for each n. 



Input level n 


3 


4 


5 


6 


7 


8 


9 


Value of c when £ = n 


20 


128 


1297 


13 660 


169 077 


2 142 197 


28 691 150 


Value of c when £ = n + 1 


8 


50 


196 


1074 


7 784 


64 528 


557 428 


Value of c when £ = n + 2 


1 


1 


1 


1 


1 


1 


1 


Final bound i?„(S'^) 


2 


2 


2 


2 


2 


2 


2 



Table 3: Running Algorithm 14.61 on the 3-sphere for all levels 2 < n < 9 



Lemma 4.7. // Algorithm \4-(>\ outputs Hn{M) < 2, then this bound is tight. In other words, 
Hn{A4) — maxi, miup excess height (p), where v ranges over all nodes at level n of and 
where p ranges over all simplification paths that begin at v. 

Proof. The result is trivial for IIn{A4) < 1. Suppose then that iJ„(A^) — 2 and this bound is not 
tight; that is, max^. miup excess height(p) < 1. When Algorithm 14.61 processes level £ ^ n + 1, it 
effectively enumerates all paths in ^i(A^) that stay between levels n and n + 1, since every arc 
between these levels emanates from one of our starting nodes at level n. Therefore Algorithm l4.6l 
would terminate at level £ < n + 1 and output iJ„(7W) < 1 instead, a contradiction. □ 
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As a result, our bounds Hn{S^) = 2 for 2 < n < 9 are all tight. Note that Lemma [4.71 does 
not work in the opposite direction: it is possible to have max^ minp excess height (p) = 2 but 
Hn{M) > 2, since the optimal path might bounce around between levels n + 1 and n + 2, using 
intermediate nodes that never appear in our subgraph G. 

It is straightforward to show that the space and time complexities of Algorithm l4.6l are linear 
and log-linear respectively in the number of nodes in G (other small polynomial factors in n and 
i also appear). Nevertheless, the memory requirements for n = 8 were found to be extremely 
large in practice (~30GB): by the time the algorithm terminated at level £ = 10, we had visited 
185 697 092 nodes in total. For n = 9 the memory requirements were too large for the algorithm 
to run (estimated at 400-500 GB), and so a two-phase approach was necessary: 

Algorithm 4.8 (Two-phase algorithm for computing Hn). The following algorithm will either 
compute the same bound Hn{A4) as Algorithm \4-(>] or will terminate with no result. Once 
again, we progressively build a finite subgraph G C J^i{A4) and track the number of connected 
components c. 

1. Initialise G to all of level n of ^i{M), as before. 

2. Construct all arcs from nodes in G at level n to (possibly new) nodes in ^i{M.) at level 
n + I. Insert these arcs and their endpoints into G, and update c accordingly. 

If c = 1 then output Hn{M) = 1 and terminate. 

3. From each node v at level n -t- 1, try all possible octahedron flips, pillow flips and prism 
flips. Let v' he the endpoint of such a flip (so v' is also a node at level n + 1). If v' ^ G 
then merge the components and decrement c if necessary. Otherwise do nothing (since v' 
would never have been constructed in the previous algorithm). 

// c = 1 then output Hn{A4) = 2 and terminate; otherwise terminate with no result. 

In essence, we use the old Algorithm 14 . 61 for the transition from level n ^ n + 1, but then we 
use Theorem 14. 5l to "simulate" the transition from level n + 1 n + 2. We only need to consider 
octahedron, pillow and prism flips in step 3 because, by Theorem 14.51 if two nodes i>,v' G G at 
level n + 1 have arcs to some common node i^" at level n + 2, then either i' and v' are related 
by such a flip, or else v and u' must already be connected in G. 

It follows that, if this two-phase approach does output a result, it will always be the same 
result as Algorithm 14. 61 The key advantage of this two-phase method is a much smaller memory 
footprint (since it does not store any nodes at level n + 2). The main disadvantage is that it 
cannot move on to level n -I- 3 if required, and so if iJ„(A^) > 2 then it cannot output any result 
at all. 

Of course by the time we reach n = 9 for the 3-sphere, there are reasons to suspect that 
Hn{S^) = 2 (following the pattern for 3 < n < 8), and so this two-phase method seems a 
reasonable — and ultimately successful — approach. For n = 9 the memory consumption was 
^^52 63, which was (just) within the capabilities of the host machine. 

4.3 Bounding path lengths 

Our next task is to compute bounds Ln{M) for a given 3- manifold M. and a given level n > 
so that, from every node at level n of ^i(AI), there is some simplification path of length 
< Ln{M.). This time we compute L„ for all n < 9 and all closed prime orientable 3-manifolds 
M. in the census, as well as for the case M. — . 
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It is infeasible to perform arbitrary breadth- first searches through £Pi{A4), and so we restrict 
such searches using two techniques: 

(i) we only search within levels n, n + 1 and n + 2, and we only explicitly store nodes in levels 
n and n + 1; 

(ii) we only visit level n + 2 when absolutely necessary, as described by Theorem 14.51 

Of course (i) is only effective if there is a simplification path of excess height < 2. The results 
of the previous section show that this is true for Ai = and n < 9, and give us hope that it 
holds for other manifolds also. Although the shortest simplification paths might pass through 
levels n + 3 or above (and so our bounds Ln{M) might not be tight), the time and space savings 
obtained by avoiding these higher levels (which grow at a super-exponential rate) are enormous. 

In the algorithm below, we refer to steps as the quantity minimised by the breadth-first 
search. Since each octahedron, pillow or prism flip involves two Pachner moves, we use a 
modified breadth-first search for weighted graphs in which some arcs are counted as one step 
and some arcs are counted as two. Such modifications are standard, and we do not go into 
details here. 

Algorithm 4.9 (Computing L„). First identify the set S of all nodes at level n of 3^i{M) 
that have an arc running down to level n — 1. Then conduct a multiple- source breadth-first 
search across levels n and n+1, beginning with S as the set of sources, where the steps in this 
breadth-first search are as follows: 

• all arcs between levels n and n + 1 of £i^i{A4), each of which is treated as one step; 

• all octahedron, pillow and prism flips at level n + 1 of ^i{M), each of which is treated as 
two steps. 

If s is the maximum number of steps required to reach any node in level n from the initial source 
set S , then output the final bound L„{M) = s + 1. If some nodes at level n are never reached, 
then output Ln{M) — oo instead. 

This time the algorithm always terminates, since levels n and n + 1 are finite. The algorithm 
is also correct, because each source node in S has a simplification path of length 1, and each 
step in the breadth-first search corresponds to a single 2-3 or 3-2 move. The space and time 
complexities are linear and log-linear respectively in the number of nodes in levels n and n + 1 
(again with further small polynomial factors in n). 

We can make some observations on the tightness of the bounds L„(7W): 

Lemma 4.10. The bound L„{A4) computed by Algorithm \4.9\ is finite if and only if every node 
at level n of ^i{A4) has a simplification path of excess height < 2. 

// this bound is finite, then it is also tight if we restrict our attention to only simplification 
paths of excess height < 2. That is, Ln{M) = max^ miuj, length (p), where v ranges over all 
nodes at level n of ^Pi{M.), and where p ranges over all simplification paths of excess height < 2 
that begin at v. 

The proof follows directly from the structure of the breadth-first search and from Theo- 
rem 14. 5[ which shows that every step up into level n -I- 2 is either part of an octahedron, pillow 
or prism fiip, or else can be replaced with a step down into level n instead. 
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For the 3-sphere, Table|3]shows how the breadth-first search progresses for each rt in the range 
2 < 71 < 9. Each cell of this table counts only nodes at level n of ^i{S^), not "intermediate 
nodes" at the higher level n + 1. For each n the algorithm outputs a final bound of i„(S''^) = 7 
or 9, proving the 3-sphere case of Theorem 14.31 



Input level n 


3 


4 


5 


6 


7 


8 


9 


Nodes in S 


3 


46 


504 


6 975 


91283 


1 300 709 


18 361866 


Nodes 2 steps from S 


8 


38 


466 


4315 


54 698 


624 144 


8 044 998 


Nodes 4 steps from S 


7 


43 


309 


2 299 


22 634 


213 345 


2 255 191 


Nodes 6 steps from S 


2 


1 


18 


71 


462 


3988 


29 054 


Nodes 8 steps from 5* 












11 


41 


Total level n nodes found 


20 


128 


1297 


13 660 


169 077 


2 142 197 


28 691 150 


Total level n nodes missing 























Final bound L„ 


7 


7 


7 


7 


7 


9 


9 



Table 4: Running Algorithm 14.91 over the 3-sphere 



Observation 4.11. The bounds Ln in Tableware tight even if we consider all simplification 
paths (not just those with excess height <2). 

Proof. For n < 7 this follows immediately from Lemma 14.101 since any simplification path of 
excess height > 3 must have length > 7. For n = 8 or 9, the bounds are likewise tight unless 
all 11 or 41 triangulations respectively that are eight steps from S have a simplification path 
of length < 8 and excess height > 3. The only way to construct such a path is using three 
2-3 moves followed by four 3-2 moves. We can enumerate all such combinations of moves by 
computer, and wc find in all 11 cases for n = 8 and in 40 of the 41 cases for rt = 9 that no such 
path exists. □ 

Table [5] shows the corresponding results for arbitrary closed prime orientable 3-manifolds. 
Each cell in this table is a sum over all such manifolds A4 with (3i{Ai) < n, except for the final 
row which lists the largest bound Ln{A4) amongst all such manifolds. In every case we have 
Ln{.M) < 17, proving the closed prime orientable case of Theorem 14.31 



Input level n 


3 


4 


5 


6 


7 


8 


9 


Nodes in S 


8 


80 


931 


11380 


151278 


2 098 537 


30 082 708 


Nodes 2 steps from S 


2 


73 


682 


7078 


80 998 


991 080 


12 579 745 


Nodes 4 steps from S 


8 


72 


483 


4163 


38 253 


391516 


4 185 195 


Nodes 6 steps from S 


1 


13 


44 


324 


2 290 


16 714 


131 545 


Nodes 8 steps from S 








12 


69 


427 


2 591 


Nodes 10 steps from S 












12 


51 


Nodes 12 steps from S 














13 


Nodes 14 steps from S 
















Nodes 16 steps from S 














1 


Total level n nodes found 


19 


238 


2140 


22 957 


272 888 


3 498 286 


46 981 849 


Total level n nodes missing 























Maximum bound L„ 


7 


7 


7 


9 


9 


11 


17 



Table 5: Running Algorithm 14.91 over arbitrary closed prime orientable 3-manifolds 
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Furthermore, by Lemma [4.101 an immediate consequence of these results is as follows. For 
any closed prime orientable 3-manifold A4, from any node at level n of where Pi{A4.) < 

n < 9, there is a simplification path with excess height < 2. This proves the remaining closed 
prime orientable case of Theorem 14.21 

It is clear from Tables S] and [5] that most nodes at level n are very close to the source set 5, 
and only a tiny fraction require all L„(A^) steps. This prompts us to consider average lengths 
of simplification paths. In particular, we consider the following two quantities: 

• CT„ — avgj, miup length (p), where v ranges over all nodes at level n of ^i{S^), and where 
p ranges over all simplification paths that begin at v; 

• 7r„ = avgj^ miup length(p) , where v ranges over all nodes at level n of all graphs 

for closed prime orientable 3-manifolds A4 with f3i{Ai) < n, and again p ranges over all 
simplification paths that begin at v. 

In other words, if we consider the shortest simplification path from each node, then ct„ is the 
average path length over all size n triangulations of the 3-sphere, and 7r„ is the average path 
length over all size n triangulations of closed prime orientable 3-manifolds. 

By aggregating the figures in Tables 2] and [SJ we can place upper bounds on (t„ and 7r„ 
respectively (since each node that is k steps from the source set S has a simplification path of 
length < A: + 1). These upper bounds are listed in Table |6] (all numbers are rounded up). 



Input level n 


3 


4 


5 


6 


7 


8 


9 


Upper bound on On 
(3-sphere triangulations) 


3.80 


2.99 


2.76 


2.34 


2.20 


2.00 


1.89 


Upper bound on 7r„ 

(closed prime orientable 3-manifolds) 


3.22 


3.16 


2.67 


2.44 


2.21 


2.05 


1.91 



Table 6: Bounding the average lengths of simplification paths 

These averages are remarkably small, and more importantly, they appear to decrease with 
n: although the worst-case paths become longer, such pathological cases become very rare very 
quickly. This has interesting implications for generic-case complexity, which we return to in 
Section El Figure [T3] gives a graphical summary of these worst-case and average-case bounds. 

4.4 Connecting minimal triangulations 

Our final algorithms examine paths that connect different nodes at the base level /3i(A^) of 
■^i(A4). Recall from Lemma [4.11 that, for most closed prime orientable 3-manifolds M., these 
are the nodes that correspond to minimal triangulations of M.. 

In particular, these algorithms compute bounds H^i^{M) and Lniin(-A^) for a given 3- 
manifold M so that, between any two nodes at level [ii{M) of ,^i(A4), there is a path of 
excess height < i7i„in(A4), and a (possibly different) path of length < Lmin(A^). 

Here the running time and memory use are less critical: because the census contains so few 
minimal triangulations of closed prime orientable 3-manifolds (just 4472 for n < 9), we can 
afford to search exhaustively through entire levels of ^i{M) as long as these levels do not grow 
too high above /3i(A^). 
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Worst case (3-sphere) 

Worst case (closed prime orientable) 

Average case (3-sphere) 

Average case (closed prime orientable) 




Level (n) 



Figure 13: Summary of bounds on lengths of simplification paths 



Algorithm 4.12 (Computing -ffmin)- before, this algorithm runs by progressively building a 
finite subgraph G C 0^i{M). We track the number of distinct components of G (denoted by c) 
and the maximum level of any node in G (denoted by i). 

1. Initialise G to all of level j3i{M) of set the number of components c to the number 
of nodes at level j3i{M), and set the maximum level £ — /3i{A4). 

2. While c > 1, expand G using a multiple- source breadth-first search as follows: 

• Maintain a queue of nodes to be processed. Initially this queue should contain all 
nodes at level £ in G (these are the multiple source nodes). 

• To process a node v, follow all arcs from v in ^i{A4) whose endpoints lie in levels 
< ^+1, and add these arcs and their endpoints to G if not already present. Each new 
node that is added to G must also be pushed onto the queue for processing. Decrement 
c each time a new arc joins two disconnected components ofG. 

• Once the queue is empty (i.e., there are no more nodes to process), increment £. 

3. Once we have c — 1, output the final bound i?min(.^) ~ ^ ^ o,''^d terminate. 

As in Algorithm 14.61 we enumerate arcs from a node u by performing 2-3 and 3-2 moves, 
and we track connected components of G using union-find. 

It is clear from the structure of the breadth- first search that, at each stage of the algorithm, 
the subgraph G contains precisely those nodes that can be reached from level (3i{M.) of ^i{A4) 
using a path that never travels above level £. We can therefore conclude: 

Lemma 4.13. The bound i?min(-^) output by Algorithm \4.12\ is tight. That is, -ffmin(-^) — 
maxj^j_^2 miup excess height(p), where and range over all nodes at level /3i(A^) of ^i{A4), 
and where p ranges over all paths in from vi to V2. 

Running this over all 1900 closed prime orientable 3-manifolds M with /3i(A^) < 9, we find 
that i?min(A1) < 2 in every case but one. The exception is one of the smallest cases in the 
census: the lens space L(3, 1), with /3i(L(3, 1)) ~ 2 and final bound i7niin(i(3, 1)) ~ 3. Table[7] 
shows a detailed breakdown of the results, which together constitute a computer proof of the 
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Base level /3i(A^) 


2 


3 


4 


5 


6 


7 


8 


9 


Total 


Manifolds with HminiM) = 


6 


7 


13 


23 


47 


106 


235 


575 


1012 


Manifolds with Hmin{M) = 1 




1 


1 


7 


20 


52 


156 


455 


692 


Manifolds with ii"min(A^) = 2 




1 




1 


7 


17 


45 


124 


195 


Manifolds with HminiM) = 3 


1 
















1 


Manifolds with HminiM) > 4 






























Table 7: Bounding excess heights of paths between nodes at level /3i{M) 



excess height results from Theorem l4.4l Note that the manifolds in the table with -ffmin(A^) ~ 
are those with only one node at level l3i{M) of ^i{M). 

To bound path lengths, we adopt a similar strategy to that for simplification paths: because 
excess heights are typically bounded above by i7,nin(-M) < 2, we perform breadth- first searches 
that stay within levels (ii{M), + 1 and l3i{M) + 2. Once again, we use octahedron, 

pillow and prism flips to avoid explicitly stepping into level /3i(A^) + 2. The full algorithm is as 
follows. 

Algorithm 4.14 (Computing Lmin)- For each node v at level Pi{Ai), run a single-source 
breadth-first search from v across levels /3i(A4) and /3i(A^) + 1 o/ 3^i{M.). The steps in this 
breadth-first search are as follows: 

• all arcs between levels j3i{M) and j3i{M) + 1, each of which is treated as one step; 

• all octahedron, pillow and prism flips at level l3i{M) + 1, each of which is treated as two 
steps. 

Let d(v) denote the maximum number of steps required to reach any node at level from 
the source v, or let d(y) = oo if some node at level /3i(A4) was never reached. After computing 
d{y) for each source node v, output the final bound Lmin(A^) — maxjy d{v). 

This approach of running a separate breadth- first search from each node v at level f3i{M) is 
perhaps wasteful, but there are so few minimal triangulations in the census that its performance 
is adequate nonetheless. 

As with Algorithm 14.91 for simplification paths, we can make some simple observations on 
the tightness of the bounds imin(A^). The following result is an immediate consequence of 
Theorem 14.51 the breadth-first search structure of Algorithm I4.14[ and the tightness of the 
height bounds iJ,„i,i(A^) as shown by Lemma [4. 131 

Lemma 4.15. The bound imin(A^) that is computed by Algorithm \4-. 14\ is finite if and only if 
HaiiniM) < 2. Moreover, if the bound LminiM) is finite, then it is also tight if we restrict our 
attention to paths of excess height < 2. That is, Lmin(A^) = T[iaxvi,v2 min^ length(p), where vi 
and V2 range over all nodes at level j3i{Ai) of !^i{M.), and where p ranges over all paths in 
^i{M) from i>i to i>2 with excess height < 2. 

Table [8] shows the results of this algorithm when run over all 1900 closed prime orientable 3- 
manifolds A4 with /3i(7W) < 9. There is just one case with L^in{M) = oo (the case A4 = L{3, 1) 
from before), and for every other manifold we have imin(A^) < 18 (note that Lmin(A^) must 
always be even). This establishes the length results of Theorem l4.4l in all cases but Ai = L(3, 1). 
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Base level Pi{M) 


2 


3 


4 


5 


6 


7 


8 


9 


Total 


Manifolds with Luiin{M) = 


6 


7 


13 


23 


47 


106 


235 


575 


1012 


Manifolds with Lmin(A^) = 2 






1 


6 


12 


33 


89 


244 


385 


Manifolds with Lmin(A^) = 4 








2 


10 


19 


53 


161 


245 


Manifolds with Luiiu{M) = 6 




1 






4 


10 


37 


88 


140 


Manifolds with Luiiu{M) = 8 










1 


5 


11 


46 


63 


Manifolds with L^in{M) = 10 












1 


7 


22 


30 


Manifolds with Lmin(X) = 12 




1 








1 


2 


7 


11 


Manifolds with Lmin{M) = 14 














1 


5 


6 


Manifolds with L^in{M) = 16 














1 


4 


5 


Manifolds with Lmin(A^) = 18 
















2 


2 


Manifolds with 18 < Lmin{M) < oo 





























Manifolds with I/min(A^) = oo 


1 























1 



Table 8: Bounding lengths of paths between nodes at level (3i{M) 



For the special case L(3, 1), there are precisely two one- vertex triangulations in the census 
with n = 2 tetrahedra; their isomorphism signatures are cMcabbgaj and cMcabbjak. A quick 
search shows that these can be connected using six 2-3 and 3-2 moves: 

2- 3 2-3 2-3 

cMcabbgaj > dLQacccbgf g > eLPkbcdddackf f > f vPQccdeedegovggo 

3- 2 3-2 3-2 

> eLPkbcdddacrkk > dLQacccbgf o > cMcabbjak 

See the appendix for details on how to "decode" these isomorphism signatures back into full 
3-manifold triangulations. Since iJmin(i(3, 1)) = 3, no fewer than six moves are possible. 

Remark. Table H] shows two "outlier" manifolds with a small base level (3i{M) = 3 but un- 
usually large bounds imin(-A4) = 6 and 12. These are the lens spaces L(4, 1) and L{5, 2), which 
are the only closed prime orientable 3- manifolds for which /3i(Ai) is not the size of a minimal 
triangulation (Lemma 14. 

The space L{4, 1) has seven one-vertex triangulations with n = 3 tetrahedra, and gives 
bounds i?min = 2 and Lmin — 12. The space L{5, 2) has five one-vertex triangulations with 
n = 3 tetrahedra, and gives bounds -ffmin = 1 and £min = 6. 

4.5 Parallelisation and performance 

For large n, the algorithms that bound simplification paths (Algorithms 14. 6[ 14.81 and 14. 9p all 
have substantial running times and memory requirements. This is due to the large number of 
nodes that they process at levels n and n -I- 1 (and in some cases, n + 2) of the corresponding 
Pachner graphs. 

Figure [U plots the actual running time and memory consumption from each of the calcula- 
tions carried out in Sections 14.21 and 14.31 To recap, these calculations were: 

• bounding excess heights by running the original Algorithm 14.61 and the two-phase Algo- 
rithm |4]8] over all 3-sphere triangulations of size n; 

• bounding path lengths by running Algorithm 14. 91 over all 3-sphere triangulations of size n, 
and over all size n triangulations of closed prime orientable 3-manifolds. 
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62 




456789 45678 
Input level (n) Input level (n) 

(a) Running times (b) Memory use 



Figure 14: Performance of algorithms from Section 2] 



Both plots use a log scale for the vertical axis. We omit n = 3 (where time was negligible) 
and do not show memory under 1 GB (where overheads dominate) . All computations were run 
on an 8-core 2.93 GHz Intel Xeon X5570 CPU with 72 GB of RAM. Algorithm gil] only shows 
n < 8, since the case n = 9 required too much memory to run (recall that this was why we 
developed the two-phase Algorithm 14.81 to replace it). The data points for Algorithm 14.91 on 
closed prime orientable 3-manifolds represent sums taken over all such manifoldsQ 

We can parallelise each of these algorithms by processing nodes simultaneously: for Algo- 
rithms 14.61 and 14.81 we simultaneously process nodes at the same level of the subgraph G, and 
for Algorithm 14.91 we simultaneously process nodes at the same distance from the source set 
iS*. However, we must be careful to serialise any updates to global structures (these include the 
subgraph G and union-find structures in Algorithms 14.61 and 14. 8[ and the queue of nodes to 
be processed in Algorithm 14. 9p . Because these global structures can grow super-exponentially 
large, we use a shared memory model on a single multi-core machine, and avoid distributed 
processing. 



Figure 14(a) shows the total CPU time summed over all threads of execution. When par- 
allelised over all eight cores, the wall time was close to 1/8 of these figures — for instance, the 
peak 46.6 days of CPU time for Algorithm 14.91 (closed prime orientable 3-manifolds, n = 9) 
represented just 5.9 days of wall time. Despite the serialisation bottlenecks, all three algorithms 
achieved over 98% CPU utilisation for the largest cases. This indicates that the task of iden- 
tifying and following arcs in the Pachner graphs (which could be parallelised) was significantly 
more expensive than querying and updating the global structures (which could not). 

For connecting minimal triangulations. Algorithms 14.121 and 14.141 are cheap in comparison: 
these required just 2.5 and 6.2 minutes of CPU time respectively to process the entire census of 
minimal triangulations, and each used under 1 GB of memory. 



^This is because Algorithm 14.91 can happily process all such manifolds together in the same run, avoiding the 
messy task of identifying beforehand which triangulations represent which specific manifolds. 
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5 Pathological cases 



In this section we first examine the most difficult triangulation to simplify from our tables, and 
we show how this pathological case corresponds to the "exceptional" brick in the Martelli- 
Petronio census |23] . Following this, we construct larger triangulations of graph manifolds 
that show how our excess height results do not generalise for arbitrary closed prime orientable 
3-manifolds. In Section [6] we return to the important case of the 3-sphere, where no such 
counterexamples are known. 

In the discussions below we use alphanumeric isomorphism signatures to identify specific 
3-manifold triangulations. The appendix of this paper includes a full specification detailing how 
such signatures encode tables of tetrahedron face gluings. Alternatively, the software package 
Regina [TD] can be used to "decode" these signatures back into 3-manifold triangulations. 

5.1 Triangulations requiring many moves 

Recall Theorem 14.31 where we show that one- vertex triangulations of closed prime orientable 3- 
manifolds of size n < 9 can be simplified using at most 17 Pachner moves. Although this bound 
is tight if we allow at most two extra tetrahedra fLemma I4.10p . the worst case is extremely 
rare: just one of the 50 778 377 triangulations in Table [5] requires all 17 moves. We denote this 
worst-case triangulation by TJnax- This triangulation is an outlier (the next- worst cases require 
just 13 moves), and we examine it here in detail. 

Tmax lias isomorphism signature jLAMLLQbcbdeghhiixxnxxjqisj , contains n = 9 tetrahedra, 
and represents a Seifert fibred space over the 2-sphere with three exceptional fibres; specifically, 
SFS (5^ : (2, 1) (3, 1) (11, —9)). It has the structure of an augmented triangular solid torus, a 
common construction for Seifert fibred spaces [U [27] : we build a 3-tetrahedron triangular prism, 
glue the triangular ends together, and attach either a layered solid torus or a Mobius band to 
each rectangular face. Figure [15] outlines the construction; see fl] for full details and notation. 



This Seifert fibred space has only one minimal triangulation, with n — 8 tetrahedra and 
isomorphism signature iLLLPQcbcgf f ghhhtsmhgosof . We denote this by T^nin- Unlike Tmax, 
the 8-tetrahedron TTnin is constructed not from prisms and layerings, but from the brick B5 
described by Martelli and Petronio in their census paper [23j . This brick is an 8-tetrahedron 
triangulation of the bounded Seifert fibred space SFS {D : (2, 1) (3, 1)), where D denotes the 
2-dimensional disc, and the three boundary edges of -B5 describe relatively complex curves on 
its torus boundary. 




Layered solid 



Mobius 
band 



Layered solid 
torus 



Figure 15: Building the augmented triangular solid torus Tn 
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The brick is unique in the Martelh-Petronio census in that it is the only large brick that 
allows Seifert fibred spaces to be built using fewer tetrahedra than standard prism-and-layering 
constructions. However, it only appears rarely in the census; this is due to the specific param- 
eters SFS {D : (2, 1) (3, 1)), the complex boundary curves, and the large number of tetrahedra. 
Indeed, amongst all minimal triangulations of closed prime orientable 3-nianifolds of size < 8, 
the triangulation %nin is the only one to contain at all. 

In a sense then, the 17 moves needed to simplify T^nax indicate that B^ is "substantially 
different" from standard prism-and-layering constructions: because T^in is the only smaller 
triangulation of the same manifold, we are forced to reorganise the tetrahedra of Tma^ to form 
a Bq configuration in order to simplify it. 

It is tempting to use B^ to search for larger triangulations that require even more moves 
to simplify, but initial attempts are unsuccessful. Stepping up to nine tetrahedra, there are 
three minimal triangulations of closed prime orientable 3-manifolds that include the brick B^: 
these represent the spaces SFS {S^ : (2, 1) (3, 1) (13, -11)), SFS {S^ : (2, 1) (3, 1) (16, -13)) and 
SFS (S*^ : (2, 1) (3, 1) (17, —14)). The corresponding augmented triangular solid tori each have 
n ~ 10 tetrahedra, and a computer search shows that each can be simplified using 17 Pachner 
moves once again. 

5.2 Triangulations requiring greater excess height 

In Theorems 14.21 and 14. 4[ we show that one-vertex triangulations of closed prime orientable 
3-manifolds of size n < 9 can be simplified using at most two extra tetrahedra, and that 
any two minimal triangulations of such a manifold can be connected using at most two extra 
tetrahedra. Here we show that for arbitrary manifolds, these results do not generalise: in 
particular, we construct triangulations of graph manifolds with ti = 10 tetrahedra for which 
three extra tetrahedra are required. 

All of the graph manifolds described in this section are obtained by joining two bounded 
Seifert fibred spaces along their torus boundaries according to a given 2x2 matching matrix. 
We refer the reader to [4] for further details and notation. 

Theorem 5.1. There exists a closed prime orientable 3-manifold M and a node v G =^i(A^) 
at level > f5i{A4) from which every simplification path has excess height > 3. 

Proof. Let M be the graph manifold SFS {D : (2, 1) (3, 1)) U/ [ -J J ] SFS {D : (3, 2) (3, 2)). The 
census of minimal triangulations in ^4j shows that A4 has just one minimal triangulation; this 
has l3i{A4) = 9 tetrahedra and isomorphism signature jLLALPQaceef gihhijkuxpwhwns. 

Let V be the node of ^i{A4) that describes the triangulation with size n = 10 and isomor- 
phism signature kLLzLQAkaceiggghij j jkxuaatlsqw. A computer search shows that no path 
from 1/ with excess height < 2 can reach level 9 of However, with excess height 3 we 

obtain a simplification path that connects v with the unique minimal triangulation described 
above. □ 

The example above was found as a side-effect of processing the 10-tetrahedron census of min- 
imal triangulations [4]: it was difficult to identify which manifold the triangulation represented, 
and this was only resolved after (with some difficulty) converting it into a known triangula- 
tion of the graph manifold A4. No other pathological cases have been found in this way. The 
10-tetrahedron triangulation above is not based on a standard prism-and-layering construction, 
although the 9-tetrahedron minimal triangulation of Ai is. 
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Theorem 5.2. There exists a closed prime orientable 3-manifold Ai and nodes t^i,J^2 G £^i{A4) 
at level /3i(A^) for which every path that joins vi with vi has excess height > 3. 

Proof. Recall that, unlike the census of all triangulations of closed prime orientable 3-nianifolds, 
the census of all minimal triangulations of such manifolds extends to 11 tetrahedra [S]. By 
running Algorithm 14.121 over this extended census, we identify seven manifolds A4 for which 
Hmin{-M) = 3 (and none for which Hj^ini-M) > 3). By Lemma [4.131 it follows that each of 
the corresponding Pachner graphs has nodes 1^1,1^2 with the property above. Table |9] 

summarises these results. □ 



Level 




Manifold M 




Nodes at 


Bound 


MM) 












level l3i{M) 




10 


SFS {D 


(2,1) (3,1)) U/ 


1 -1 
1 


SFS {D : 


(3,2) (3,2)) 


6 


3 


11 


SFS(D 


(2,1) (3,1)) U/ 


-1 
1 - 


] SFS(D 


:(3,2) (3,2)) 


18 


3 


11 


SFS(L> 


(2,1) (3,1)) U/ 


1 -1 
1 


SFS {D : 


(3,2) (4,1)) 


9 
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11 


SFS(L> 


(2,1) (3,1)) U/ 


1 -1 
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(3,2) (4,3)) 
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11 


SFS (D 


(2,1) (3,1)) U/ 


1 -1 
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SFS (D : 


(3,2) (5,2)) 


9 
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11 


SFS (D 


(2,1) (3,1)) U/ 


1 -1 
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SFS (D : 


(3,2) (5,3)) 


9 


3 


11 


SFS (D 


(2,1) (4,1)) U/ 


1 -1 

1 


SFS {D : 


(3,2) (3,2)) 


6 
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Table 9: Manifolds from the 11-tetrahedron census with H^in{M) > 2 



As with Theorem 15.11 all of the manifolds in Table [9] are graph manifolds. We return to 
this issue in Section [51 where we discuss the lack of known pathological examples for "simple" 
manifolds (in particular, the 3-spliere). 

6 Discussion 

As noted already, the bounds on simplification paths that we obtain in Section[3|are astonishingly 
small. Although we only consider triangulations of size n < 9, this is not a small sample: 
the census includes ^ 150 million triangulations, including ^ 31 million one-vertex 3-spheres, 
plus another ^51 million one- vertex triangulations of 1 900 distinct closed prime orientable 
3-manifolds that cover all eight Thurston geometries as well as non-geometric combinations of 
these [24]. 

These results have interesting implications for simplification algorithms. Given a triangula- 
tion T (of arbitrary size) that wc wish to simplify, we can follow the strategy of Algorithm 14.91 
and perform a breadth-first search from T using Pachner moves and octahedron, pillow and 
prism flips, never adding more than two extra tetrahedra. Theorem 14.21 gives us hope that this 
will indeed simplify T, and Theorem 14.31 gives us hope that the search will be relatively fast. 
Initial observations from using such simplification techniques in ongoing research projects show 
them to be remarkably successful. 

Triangulations of the 3-sphere exhibit particularly good behaviour. In contrast to closed 
prime orientable manifolds, the worst-case number of moves needed to simplify a 3-sphere trian- 
gulation barely grows at all with the size n (Figure [T3|. Moreover, unlike the graph manifolds 
discussed in Section [5l there are no known pathological 3-sphere triangulations that require 
more than two extra tetrahedra to simplify, despite the 3-sphere having far more triangulations 
than any other manifold in the census. These results lead us to make the following conjecture: 
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Conjecture 6.1. From any node at any level n > 3 of the graph 3^i{S^), there is a simplification 
path of excess height < 2. 

If true, this would help explain why 3-sphere triangulations are so easy to simplify in practice. 
Moreover, by combining this result with Theorem 13. 1[ we could reduce Mijatovic's bound on 
the number of Pachner moves needed to simplify a one- vertex triangulation of the 3-sphere from 
exp(0(n^)) down to exp(0(ri log n)) instead]! 

There are theoretical reasons to believe that Conjecture 16.11 might be true. Similar results 
are known for the simplification of knots in R'^: Dynnikov shows that an arc presentation of 
the trivial knot can be simplified using elementary moves without increasing the complexity of 
the diagram |131 114| . and Henrich and Kauffman use his results to show that any projection of 
the unknot can be simplified using Reidemeister moves with only a quadratic increase in the 
number of crossings [16]. In the setting of triangulation simplification, Mijatovic's bounds for the 
3-sphere — though extremely large — remain far smaller than for other manifolds (whose bounds 
involve towers of exponentials). Moreover, every pathological triangulation from Section [5] that 
requires three extra tetrahedra is of a non-geometric graph manifold: it is possible that the 
boundaries between geometric components in such manifolds act as barriers to simplification 
that would not exist for the 3-sphere. 

We now turn our attention to path length; that is, the number of Pachner moves required to 
simplify a triangulation. Our next conjecture involves the average simplification path lengths 
an (for 3-sphere triangulations) and 7r„ (for closed prime orientable 3- manifold triangulations), 
as defined in Section lOl 

Conjecture 6.2. Both quantities cr„ and 7r„ are in 0(1); that is, they are bounded above by a 
constant that does not depend on n. 

This conjecture is clearly supported by the results in Figure [131 Looking more closely at 
Tables m and [5] however, we can formulate a more interesting observation: 

Conjecture 6.3. For a given 3-manifold M., let ipniM.) denote the fraction of nodes at level n 
of S^i{AA) that have an arc leading directly down to level n ~ 1. Then, if M is the 3-sphere or 
a closed prime orientable 3-manifold, lim„_^oo 0n(-^) = 1- 

In other words, as n ^ cxo, almost all size n triangulations of the 3-sphere or of closed prime 
orientable 3-manifolds can be simplified immediately, using just a single 3-2 move. Intuitively 
this seems reasonable, since with more tetrahedra it is more likely that a suitable degree three 
edge can be found. However, proving it appears difficult: like most probabilistic arguments 
involving 3-manifold triangulations, we run into the critical problem that almost all pairwise 
gluings of tetrahedron faces do not represent 3-manifolds at all [12] . 

This observation can be framed in terms of generic complexity [33], where we are allowed 
to ignore a vanishingly small population of pathological inputs. In this context. Conjecture 16.31 
states that generic triangulations of or of closed prime orientable 3-manifolds can be simplified 
in a single move. An interesting direction for future research is to build a 3-sphere recognition 
algorithm based on Pachner moves that runs in polynomial time for generic inputs. 

We return now to worst case path length for 3-sphere triangulations. Again, we observe 
from Figure [TSl that the worst-case bound Ln{S^) barely changes at all for 3 < n < 9: there 
is just one increment, which is of the smallest possible size (recalling that L„ must always be 

*Even if the worst-case excess height for simplification paths in grows at a rate of 0(n), we can still 

reduce Mijatovic's bound to exp(0(n log n)) in this way. 
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odd). This is in contrast to the case of arbitrary closed prime orientable manifolds, where the 
worst-case path length grows with n, and does so at a clearly accelerating rate. 

This slow rate of growth for Ln{S^) has direct implications for the complexity of 3-sphere 
recognition: 

Observation 6.4. If Ln{S^) € o{n/ \ogn), then there is a 3-sphere recognition algorithm with 
running time in o(a") for any a > 1. That is, the 3-sphere can be recognised in sub- exponential 
time. 

Proof. From any n-tetrahedron triangulation T there are at most 0{n) possible 2-3 or 3-2 moves 
(one for each face or edge of T) ■ We can therefore enumerate all possible sequences of k moves 
in O (n* • {n -f- fc)*^) time, where n* is a small polynomial term that accounts for the mechanics 
of performing each 2-3 or 3-2 move. Setting k = Ln{S^) and observing that n + k e 0{n), we 

find that if T can be simplified, this can be done in O ^n*"'"^"^'^^-'^ time. 

Repeating this operation n — 2 times, we can either simplify 7" to a known 2-tetrahedron 
triangulation of S^, or else show this to be impossible. The total running time is 

O (n ■ n*+^"('5')) = O (n*+i • exp (logn ■ i„(S'^))) . 

Since L„(S'^) €E o(n/ log n), this running time is o(a") for any a> 1. □ 

Although there is too little data to predict the precise growth rate of Ln{S^) (a situation 
that Theorem 13 . 1 1 shows is extremely difficult to rectify), it is certainly plausible to suggest from 
Figure fT3l that Ln{S^) is "sufficiently sublinear" to satisfy Ln{S^) & o{n/\ogn). This leads us 
to the following conjecture: 

Conjecture 6.5. From any node at any level n > 3 of the graph 3^i{S^), there is a simplification 
path of length < Ln{S^) where Ln{S'^) G o(n/ log n), and therefore sub- exponential time 3-sphere 
recognition is possible. 

If this conjecture could be proven, the resulting sub-exponential time 3-sphere recognition 
algorithm would be a significant breakthrough in algorithmic 3-manifold topology. 
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Appendix: Isomorphism signatures 

In Section [3.21 we describe isomorphism signatures, which are small pieces of data that uniquely 
define the isomorphism class of a triangulation, and that can be computed in small polynomial 
time. Recall that to compute an isomorphism signature, we must: 

(i) enumerate all 0{n) canonical labellings of the input triangulation; 
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(ii) encode the full set of face gluings for each canonical labelling; 

(iii) choose the lexicographically smallest encoding amongst all canonical labellings. 

In this appendix we specify the encoding used in step (ii). For the theoretical arguments in 
Section 13.21 it is sufficient to use a simple bit sequence that encodes the complete table of face 
gluings. In practice however, we use a more compact method: 

• we do not encode redundant information from the face gluings table; 

• we use printable characters (letters, digits and punctuation) instead of plain bits. 

In this way, the encoding is kept short but remains easy to write or type by hand. This is useful 
in both papers and software: authors can describe triangulations precisely without presenting 
large tables of face gluings (as seen in Section [5] of this paper), and readers can reconstruct these 
triangulations quickly using the appropriate software. 

The precise encoding described here is implemented in version 4.90 of the freely-available 
software package Regina [lOj . This encoding method is defined not only for closed triangulations 
(as described in this paper), but also for triangulations with boundary (where some tetrahedron 
faces are not glued to any partner face) and ideal triangulations (where some vertex neigh- 
bourhoods are not bounded by spheres). The encoding draws upon and extends ideas from the 
dehydration format of Callahan et al. [TT] . 

Encoding integers 

A key piece of the encoding process is to translate integers into printable characters. We describe 
two methods: one for "small" integers in the range 0, . . . , 63, and one for "large" integers in the 
range 0, . . . , n. 

For small integers, we use the following translation table from 0, . . . , 63 to printable char- 
acters. This is similar to, but not the same as, a typical Base64 encoding tablelfl For each 
i € {0, . . . , 63} we let 7r(z) denote the corresponding printable character. 



Integer i 


•• 


• 25 


26 •• 


• 51 


52 •• 


• 61 


62 


63 


Printable character 7r(i) 


a • • 


z 


A • ■ 


Z 


• • 


• 9 


+ 





To encode large integers we use d printable characters, where d = [logg4(n)J -1-1. Specifically: 
for any integer i G {0, . . . ,n}, we write i as a d-digit number in base 64, and then use 7r(-) to 
encode the individual digits as described above, beginning with the least significant digit. We 
denote the resulting string of characters by e(i). 

For example, suppose that n = 100 and we wish to encode the integer i — 93. This encoding 
uses d = [log64(100)J +1 = 2 characters. Since 93 = 29 4- 1 • 64, the two base 64 "digits" of i 
are (29, 1). The final encoding is then e(93) = 7r(29)7r(l) = Db. 

If j n, we may need to pad the base 64 representation of i with leading zeroes. For 
instance, if n = 100 and we wish to encode i — 5, the two base 64 "digits" are (5,0), and the 
final encoding is e(5) = 7r(5)7r(0) = fa. 

key difference is that we use lower-case letters first, to simplify the process of manually typing isomorphism 
signatures on the keyboard. 
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Encoding face gluings 

Here we identify the information required to reconstruct the full table of face gluings for a 
canonical labelling. Such a table can be seen in Table [2] (from Section |3^ . which we use as an 
example throughout this discussion. 

We begin by removing redundant information from the gluings table. Suppose that face / 
of tetrahedron t is glued to face /' of tetrahedron t' . If (t', /') < {t, f) lexicographicallj0, then 
we delete the table cell in position {t,f). This is because the same gluing has already been 
seen from the other direction in cell {t',f'), and so the information in cell {t,f) is redundant. 
Table [TOl shows a copy of Table [2] with these redundant cells crossed out. 





Face 
Vertices 123 


Face 1 
Vertices 023 


Face 2 
Vertices 013 


Face 3 
Vertices 012 


Tet. 
Tet. 1 
Tet. 2 


Tet. 
Tet. 2 
Tet. 2 


120 
230 
012 


Tet. 1: 023 


Tet. 2: 013 
Tet. 1: 012 


TM.-er-3T2 
Tit,-4r-0Tf 



Tabic 10: Removing redundant information from the face gluings table 



To reconstruct the table cells that remain, we extract the following information: 

• Destination sequence: As in Section [3. 21 let At.f denote the tetrahedron glued to face / of 
tetrahedron t, so that Atj G {d, 0, . . . , n — 1} for alH = 0, . . . , n — 1 and / = 0, . . . , 3. Here 
we use the special symbol Atj = d to indicate that face / of tetrahedron t is a boundary 
face, i.e., it is not glued to any partner tetrahedron face at all. The destination sequence 
is obtained by writing out ^o,Oj ^o,i, ^o,2j ^o,3j ^i,0: • ■ • i ^n-i,3i and then deleting any 
entries Atj that correspond to deleted table cells (<,/). 

For our example in Table [TUl the corresponding destination sequence is 0, 1,2, 2, 1,2. 

• Type sequence: For each term At j in the destination sequence, we classify the correspond- 
ing tetrahedron face as one of the following three types: 

— Type 0: This indicates that the face is a boundary face, i.e., Atj ~ d. 

— Type 1 : This indicates that the face is joined to a new tetrahedron. That is, Atj = k 
where (i) k > 1, and (ii) the term Atj is the first time that k appears in the 
destination sequence. 

— Type 2: This indicates that the face is joined to an old tetrahedron. That is, either 
Atj = 0, or else Atj = k > 1 but there is some other term Af.f = k that appears 
in the destination sequence earlier than Atj. 

We obtain the type sequence by writing out the corresponding face type for each term in 
the destination sequence. 

For our example, the type sequence is 2,1,1,2,2,2. The two 'I's in this sequence corre- 
spond to the first appearances of tetrahedra 1 and 2, in cells (0, 1) and (0, 2) respectively 
of the face gluings table. 

'That is, either t' < t, or else t' = t and /' < /. 
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• Permutation sequence: For each non-boundary tetrahedron face {t, f) that is glued to some 

partner face {t',f'), wc identify the corresponding gluing perm,utation ptj € S^. This is 
the permutation for which Ptj{f ) = /', and each vertex v ^ f oi tetrahedron t maps to 
vertex Ptj{v) ^ f of tetrahedron t' under this face gluing. 

We can represent each permutation in ^4 by a number: simply order them lexicograph- 
ically and number them 0, . . . , 23. The identity (0,1,2,3) M> (0,1,2,3) becomes 0, the 
permutation (0,1,2,3) h-t- (0,1,3,2) becomes 1, and so on until (0,1,2,3) ^ (3,2,1,0) 
becomes 23. 

To build the permutation sequence^ we write out the number for each permutation p^j 
corresponding to each term At.f ^ d the destination sequence. 

In our example, this sequence is 21, 0, 0, 9, 1, 18. The two 'O's come from the faces of type 1: 
because the labelling is canonical, those permutations must be the identity. 

The full isomorphism signature 

To completely encode a labelled triangulation, we bundle all of the above information into a 

sequence of printable characters as follows. 

• We begin the sequence by encoding n and d (recall that d indicates how many printable 
characters are required to encode each "large integer" in the range 0, . . . ,n). Specifically: 

if n > 63 then we begin with the special marker 7r(63) followed by 7r((i) and e(n). If n < 63, 
we simply begin with 7r(n) and it is implicitly understood that d = \. 

Strictly speaking, because we encode 7r((i), this scheme only works \id < 64. However, this 
is true for all n < 64^^ ~ 6 x 10^^^, which covers all conceivable practical requirements. 

• We follow this with the type sequence. Since each face type is in {0, 1,2}, we can encode 

three terms of the sequence in each printable character. Specifically, if the type sequence 
is to, ti,. . ., then we write 7r(to + 4ii + 16^2), 7r(i3 + 4^4 + 16^5), . . . , padding the sequence 
with one or two trailing zeroes if necessary. 

• Next we write the destination sequence, but only for faces of type 2. Each destination 
At J is encoded as €{Atj). We ignore the other faces because type faces are boundary, 
and type 1 destinations can be deduced from the fact that the labelling is canonical. 

• We finish with the permutation sequence, again only for faces of type 2. Each permutation 
number i is encoded as 7r(z). We ignore type 1 faces because, in a canonical labelling, all 
type 1 gluings must use the identity permutation. 

For our example: since n = 3 < 63, we begin the sequence with 7r(3) = d. The type sequence 
2, 1,1, 2, 2, 2 encodes to 7r(2 + l-4 + l-16) followed by 7r(2 + 2-4 + 2-16); that is, 7r(22)7r(42) = wQ. 
The destination sequence for type 2 faces is 0, 2, 1, 2, which encodes to acbc. The permutation 
sequence for type 2 faces is 21,9,1,18, which encodes to vjbs. We combine these pieces to 
obtain the final sequence dwQacbcvjbs. 

To obtain an isomorphism signature, we must choose the lexicographically smallest encoding 
amongst all canonical labellings. For this we order printable characters according to the standard 
ASCII tables (so, for instance, 1 < Z < a). This keeps things simple for programmers, who can 
using ready-made comparison functions such as strcmp in C. 
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In our example, the encoding dwQacbcvjbs from Table [TUl is not an isomorphism signature. 
Several different canonical labellings give lexicographically smaller encodings; the smallest is 
dLQabccbcjj, which becomes the isomorphism signature for this triangulation. 

As a final remark, we can easily extend isomorphism signatures to support disconnected 
triangulations (which are not considered in this paper): simply list the signatures for each 
connected component in sorted order. 
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